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Abstract 

In this paper we present a measurement of the natural widths of and 
Using data from the FOCUS experiment, we find r(S°) = 1.55+[j ± 0.38 MeV/c^ 
and r(S++) = 2.05+o;38 ± 0.38 MeV/c^. The first errors are statistical, the second 
systematic. These results are obtained with a sample of 913 — > A^vr" decays 
and 1110 A^vr^ decays. These results are compared with recent theoretical 

predictions. 

PACS numbers: 14.20.Lq 13.30Eg 



The natural widths of the excited charmed baryons have become experimen- 
tally accessible [1-3] only in the last several years. The widths of most of the 
known excited charmed baryons are poorly determined with either no exist- 
ing measurement or with an upper limit ;Q most of the lower mass excited 
states appear to have widths comparable to or less than the resolution of cur- 
rent experiments. Combined with relatively low statistics, this has made such 
measurements challenging. Direct measurements of the T,c widths are impor- 
tant since most of the current theoretical models predict the widths of charmed 
baryons by extrapolating from the hyperon widths. Accurate measurements 
of the Sc widths will enable more accurate predictions of the widths of these 
and other excited charm states and test the physics underlying these models. 

In this paper, we use data from the FOCUS experiment to obtain measure- 

^ During final preparation of this manuscript, we became aware of a measurement 
of and S^"*" intrinsic widths by the CLEO collaboration [4]. 
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ments of r(S°) and r(S^"'"). The and mass differences with respect to 
the were presented in an earher paper [5]. FOCUS is an upgraded version 
of FNAL-E687 [6] which collected data using the Wideband photon beam- 
line during the 1996-1997 Fermilab fixed-target run. The FOCUS experiment 
utilizes a forward multiparticle spectrometer to study charmed particles pro- 
duced by the interaction of high energy photons {{E) ^ 180 GeV) with a 
segmented BeO target. 

Charged particles are tracked within the spectrometer by two silicon microver- 
tex detector systems. One system is interleaved with the target segments; the 
other is downstream of the target region. These detectors provide excellent 
separation of the production and decay vertices. Further downstream, charged 
particles are tracked and momentum analyzed by a system of five multiwire 
proportional chambers and two dipole magnets with opposite polarity. Three 
multicell threshold Cerenkov detectors are used to discriminate among elec- 
trons, pions, kaons, and protons [7]. 

The Sc candidates are reconstructed via the decay chain Sc A'^tt'^, with 
A+ candidates found in four decay modes|^ pK'-K^ ^ pK% A°7r"*", and 

A°7r'^'7r~7r"'". The reconstruction of K^s^s and A'^'s is described elsewhere [8]. Due 
to topological differences and varying levels of background, the values of the 
analysis cuts for each of the four A+ decay modes vary. A decay, or secondary, 
vertex is formed from selected reconstructed tracks; the momentum vector of 
this charm candidate is then used as a seed to intersect other reconstructed 
tracks in the event to find the production, or primary, vertex [6]. We impose 
a minimum detachment cut which requires that the measured separation of 
these two vertices divided by the error on that measurement be greater than 
our cut (typically 3-5). We also ensure that both vertices are well formed 
by requiring a confidence level greater than 1% on the fit to each vertex. 
To remove longer lived charm backgrounds we also require the reconstructed 
proper lifetime to be less than some amount, typically 4-5 times the mean 
A^ lifetime. Finally, momentum and particle identification cuts are applied to 
each decay mode. Invariant mass plots for each of the decay modes are shown 
in Fig. 1. 

The Sc candidates are reconstructed by combining the A^ candidates within 
approximately 2 a of the mean A+ mass with a charged pion.[^ The vertex 
formed by the A+ candidate, the pion candidate, and at least one other track 
must have a confidence level greater than 1%. Cerenkov identification on the 
soft pion requires that the pion hypothesis is not heavily disfavored with re- 
spect to any other identification hypothesis. 

To remove systematic effects due to the reconstruction of the A+ mass, we 
^ Throughout, charge conjugate states are implied. 

^ Referred to as a "soft pion" since it is usually a low momentum particle. 
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Fig. 1. Invariant mass distributions for candidates. Shown are the four decay 
modes used in the analysis and, on the bottom, the combined invariant mass plot 
for all four decay modes. 

compute and plot the invariant mass differences (AM = M(A+7r'^) — M(A^)). 
Because most of the uncertainty in the mass difference arises due to multi- 
ple scattering of the soft pion, we improve the measurement of the soft pion 
momentum as follows. The primary vertex is refit without the soft pion and 
the soft pion is constrained to originate from this new primary vertex. If the 
confidence level of this constraint is less than 1%, the candidate is discarded. 
The measured A+ momentum and mass are combined with the constrained 
pion momentum and known mass to form M(A+7r^). The computed A+ mass 
is subtracted to obtain the invariant mass difference, AM. The reconstructed 
mass difference distributions are shown in Fig. 2. 

To measure r(Ec), we fit the Ec signal distributions with a relativistic constant 
width Breit-Wigner function convoluted with a parameterization of the exper- 
imental resolution derived from Monte Carlo. The experimental resolution is 
determined by generating Monte Carlo events with r(Sc) = and fitting the 
resulting Sc distributions to a sum of two Gaussians, Yield, cr, x), with 
the same central value x: 



G (AM; (1 - (71, Mfit) + G (AM; fY, a-,, Mfit) 



(1) 
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Fig. 2. Invariant mass differences for (left) and (right) candidates. The fit 
functions, over their range of 140-200 MeV/c^, are shown as are the widths and 
statistical errors. 
Table 1 

Sc resolution parameters calculated in Monte Carlo. 







0-2 


/ 


State 


(MeV/c2) 


(MeV/c2) 


Yields / (Yieldi +Yield2) 




0.84 ± 0.03 


1.93 ±0.05 


0.53 ± 0.03 




0.83 ± 0.03 


1.92 ± 0.05 


0.56 ± 0.03 



We use a background function described by 

A^(l + a{AM - m^)AM'^) (2) 

where m^r is the tt^ mass. N, a, and /3 are allowed to vary. Three parameters 
of interest are extracted from the fit: ai, (T2, and / = Y2/Y which describe, 
respectively, the resolution of the narrow and wide portions of the resolution 
function and the ratio of the second yield to the total. The extracted param- 
eters are shown in Table 1. Mat is the Ec — A+ mass difference. 

We determine the shape parameters [a and [3) of the background function 
(Eq. (2)) by fitting the sum of the M(A+7r+) and the M(A+7r-) data distri- 
butions to Eq. (2) and a convoluted Breit-Wigner signal shape. The values of 
a and /3 are fixed to these values when fitting the individual distributions. 

While the fits for and distributions to a Breit-Wigner convoluted 

with the resolution function use a total of nine parameters, only four of these 
are free parameters (yield, M(Ec — A+), r(Ec), and N, the background nor- 
malization). The parameters a, cti, (T2, and / are previously determined 
and arc fixed. Performing the fit as described above we obtain values of 
r(EO) = l.bbt^ofj MeV/c2 and r(E++) = 2.05^[J;|^ MeV/c^ with yields of 
913 ± 77 and 1110 ± 83 respectively. The extracted mass differences are con- 
sistent with our previous measurement [5]. 
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Table 2 

Comparison of D*~^ line shape parameters for data and Monte Carlo. Statistical 
errors only. 







0-2 




State 


(MeV/c2) 


(McV/c2) 


/ 


D*+ Data 


0.66 ± 0.00 


3.5 ±0.2 


0.05 ± 0.02 


D*+ MC 


0.61 ±0.00 


2.8 ±0.5 


0.05 ± 0.02 



The largest systematic uncertainties in this measurement arise from our im- 
perfect knowledge of the experimental mass resolution. To better understand 
our experimental resolution in a kincmatically similar situation, we have used 
the decay D*^ D°7r+ as a benchmark. The D*"*" is much closer to decay 
threshold and should expose any problems with our simulation of the experi- 
mental resolution. We compare the resolution of classes of events seen in the 
data with that seen in the simulation. The Monte Carlo uses an input D*'^ 
natural width of 96 keV/c^, the value from a recent CLEO [9,10] report. For 
both Monte Carlo and data we fit the distributions to the sum of Eq. (1) and 
Eq. (2). No value of r{D*~^) is determined by or included in the fit. (Since 
r(D*+) < a(AM(D*+)), the precise value of r(D*+) influences the fit results 
only shghtly. The Breit-Wigner modifications to the line shape are negligible 
and arc safely ignored.) By studying any discrepancies in the D*'^ case, we 
understand our experimental resolution better for the Ec. 

We find that for the D*"^, our overall experimental resolution matches that 
predicted by the Monte Carlo quite well and find no definitive evidence for 
any discrepancy except for cxi which is about 10% less in Monte Carlo. The 
values of the fits to the D*^ line shape for data and Monte Carlo are shown in 
Table 2. To be conservative, we test the effects of varying the experimental res- 
olution by the maximum allowed by our D*'^ studies. In our description of the 
predicted resolution we independently vary the three parameters ai, a2, and 
/ by ±10%, ±20%, and ±15% respectively. These changes induce variations 
of approximately =f0.20, =f0.25, and =f0.20 MeV/c^ in r(Ec), respectively. 

In addition, we have studied the resolution function for the D*'^ as a function 
of target configurations, the production target segment, decay modes, and the 
detachment cut. The resolution is clearly dependent on the production target 
segment and the several target configurations used in FOCUS. The dependence 
of the resolution on these two variables is well predicted by the Monte Carlo 
simulation and the distribution of Eg events as a function of these variables is 
also well matched by the Monte Carlo. The resolution function also depends 
somewhat on the kinematic variables, most strongly on the momentum of the 
excited state. As an example, in Table 3, we show the variation of o"i for 
the D*~^ in four bins of momentum. The data and Monte Carlo disagree by 
less than 10% as shown in Table 2, but the disagreement is constant vs. D*'^ 
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Table 3 

Comparison of D*~^ core resolution, o"i (MeV/c^), vs. momentum. 



p(D*+) GeV/c 


< 60 


60-80 


80-100 


> 100 


D*+ Data 
D*+ MC 


0.61 ±0.01 
0.57 ±0.00 


0.63 ± 0.01 
0.59 ± 0.00 


0.65 ±0.01 
0.62 ± 0.00 


0.70 ± 0.01 
0.66 ± 0.00 




150 200 
GeV/c 



Fig. 3. Momentum spectrum of S|? and 2++ (summed) for data (points) and Monte 
Carlo (dashed curve). Both Monte Carlo and data are sideband subtracted. 

momentum. In Fig. 3, we show that the momentum distribution is well 
modeled by the Monte Carlo. Other variations in resolution are found as a 
function of soft pion momentum and charm and soft pion directions, with a 
weaker dependence on these kinematic variables. 

We also perform two additional tests of our fitting and reconstruction method. 
In the first, we produce a large number of distributions which are statistically 
similar to the data. These distributions are then fit with the same methods 
used for the real data. The extracted fit parameters are compared with the 
known input parameters. From this source we estimate a maximum bias in 
the fitting technique of 0.05 MeV/c^ and confirm the validity of our statistical 
errors. 

In the second test, we generate Monte Carlo events with a natural width which 
is comparable to what is measured and measure the width of these events using 
the same methods. From this study we find no evidence for a systematic error. 

Adding all the systematic errors in quadrature, we find a total systematic 
error of ±0.38 MeV/c^ for each r(Ec) measurement. This gives final values of 
r(SO) = l.55l^ofj ± 0.38 MeV/c^ and r(E++) = 2.05tg:^^ ± 0.38 MeV/c^. 

There are several recent theoretical predictions for the widths of the states 
[11-16]; all predict r(Sc) in the range 1-3 MeV/c^. Several different theoreti- 
cal models are used, including the Relativistic Three Quark Model (RTQM), 
Heavy Hadron Chiral Perturbation Theory (HHCPT), and the Light Front 
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Table 4 

Predicted Sc widths from various theoretical models. All units are MeV/c^. 



Author 


r(sO) 


r(s++) 


Method 


Ivanov, et al. [11,12] 


2.65 ±0.19 


2.85 ±0.19 


RTQM 


Tawfig, et al. [13] 


1.57 


1.64 


LFQM 


Huang, et al. [14] 


2.4 


2.5 


HHCPT 


Pirjol, et al. [15] 


1.0-3.0 


1.1-3.1 


HHCPT 


Rosner [16] 


1.32 ± 0.04 


1.32 ± 0.04 


Ratios 



Quark Model (LFQM). 

These models predict partial widths, but since Eg A+tt is the only allowed 
strong or electromagnetic decay mode for the states under study, wc safely take 
the partial width as the total width for each state. The resulting predictions 
from the models are shown in Table 4. Our measurements are not of sufficient 
precision to strongly favor any one of these models over another. 

In conclusion, we present measurements of the natural widths of the S° and 
excited charmed baryons using data from the FOCUS experiment. We 
find r(EO) = l.55toi\ ± 0.38 MeV/c^ and r(S++) = 2.051^;^^ ± 0.38 MeV/c^ 
which are consistent with most current theoretical predictions. 
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collaborating institutions. This research was supported in part by the U. S. 
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Cientifico e Tecnologico, CONACyT-Mexico, the Korean Ministry of Educa- 
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